Can topological defects be formed during preheating ? 
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We study the dynamics of a scalar field $ with the potential <?(|$| 2 — n 2 ) 2 /2 (g = self-coupling 
constant, n — symmetry breaking scale) after inflation and make clear whether topological de- 
fects can ever be formed during preheating. In particular, we pay attention to GUT defects 
(77 ~ 10 15 GeV — 10 17 GeV), and consider three types of fluctuations. The first one is produced 
due to parametric resonance, the second is due to the negative curvature of the potential, and the 
last is created during inflation. We search for the parameter region that nonthermal fluctuations of 
the scalar field produced through the parametric resonant decay of its homogeneous part do not lead 
to defect formation. We find that this region is rather wide, and the GUT defects are not produced 
after inflation. This fact shows that the positiveness of the effective mass square of the field and 
production of large fluctuations whose amplitude is as large as that of homogeneous mode are not 
enough conditions for full symmetry restoration. 



PACS numbers: 98.80.Cq, latexll.27.+d 



I. INTRODUCTION 

The inflation models Glp | were invented to solve several 
cosmological problems such as the flatness problem and 
horizon problem. Among many inflation models it seems 
that the most natural and simplest model is chaotic in- 
flation in which one only needs a flat potential for 
a scalar field called the inflaton. The quantum fluctua- 
tions of the inflaton field during inflationary epoch be- 
come the density fluctuations with scale-invariant spec- 
trum which accounts for the structure formation of the 
universe. The amplitudes of the fluctuations are deter- 
mined by the self-coupling constant A of the inflaton field 
6 (for the V = (A/4> 4 model), which should be ~ 1CT 13 
j| to explain the anisotropics of Cosmic Microwave Back- 
ground (CMB) radiation observed by the Cosmic Back- 
ground Explorer (COBE) §. 

In order that the radiation-dominated universe should 
be recovered after inflation, the reheating process is nec- 
essary for transferring the vacuum energy to relativistic 
particles. The inflaton field that is oscillating at the bot- 
tom of its effective potential decays into other lighter 
particles due to its coupling to other particles. In the 
old version of reheating theory H , the decay rate is es- 
timated by using perturbation theory (Born approxima- 
tion), and the reheating temperature can be estimated 
as Trh — 10 _1 yTtot-Mp 0, based on the single parti- 
cle decay. Here r to t is the total decay rate and M p the 
Planck mass. 

Recent investigations revealed that the explosive de- 
cay of the inflaton field takes place in the first stage of 
reheating (8 17 due to the effect of parametric resonance. 
There are three stages in the reheating process: Inflaton 
field decays drastically into bose particles with broad res- 
onance. This first stage is called preheating || . After the 
back reactions of produced particles become significant, 
the resonant decay goes into the narrower band. At the 
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last stage the scatterings and further decays of created 
particles occur and the thermal equilibrium is achieved, 
which is the completion of the reheating process. 

Particle production can be described by a Mathieu- 
type equation or Lame-type equation, which has an un- 
stable solution in some regions of parameters (instability 
bands). If the relevant parameters stay in the instabil- 
ity bands long enough, the solution grows exponentially. 
This means that the number of particles produced be- 
comes exponentially large so that the parametric reso- 
nant decay of the inflaton takes place very efficiently. 
The physical meaning of this phenomenon is considered 
as an induced effect, in the sense that the presence of pro- 
duced particles stimulates the inflaton field to decay into 
those particles [0. This phenomenon is thus peculiar 
to those particles that obey Bose- Einstein statistics, i.e., 
bosons. This is the reason that the parametric resonant 
decay into fermions cannot occur due to Pauli's exclusion 
principle. 

Bosons produced at the preheating stage are far from 
thermal equilibrium and have enormously large occupa- 
tion numbers in the low-energy side of their spectrum. 
One of surprising results of these nonthermal fluctuations 
(created particles) during preheating is the symmetry 
restoration which may lead to the formation of topologi- 
cal defects |L8j (see also p9] ) . Assuming that the inflaton 
potential is given by (A/4j(0 2 — rj 2 ) 2 , the inflaton <j> decays 
into cj> particles through parametric resonance and pro- 
duced particles result in large fluctuations (S<fi 2 ) ^> n 2 , 
the effective potential of the inflaton is changed such as 
it has no negative mass, i.e., the symmetry is restored. 
Then, the restored symmetry is spontaneously broken 
again to produce topological defects. If the symmetry- 
breaking scale is the grand unified theory (GUT) scale, 
GUT defects are produced. The GUT monopoles or do- 
main walls created after inflation destroy the advantage 
of introducing the inflationary universe. On the other 
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hand, the GUT cosmic string may be another seed for 
the large scale structure formation. In Ref. fisj] , however, 
the details of the dynamics of the inflaton field were not 
studied. 

In this paper, we study the dynamics of both real and 
complex scalar field which has a flat potential V(Q) — 
|(|<I>| 2 — rj 2 ) 2 and is responsible for topological defects, 
using numerical calculations including the effect of the 
resonant decay and taking into account the back reac- 
tions of fluctuations during preheating in an expanding 
universe. Naively, the topological defects are formed for 
this type of potential if there are very large fluctuations of 
$ (| <5$ | ^S> \&\) which implies the full symmetry restora- 
tion. As mentioned above, in Ref. it is pointed out 
that these fluctuations are produced by the effect of para- 
metric resonance. Besides, it will be seen later that the 
fluctuations also grow efficiently by the effect of the neg- 
ative curvature of the potential of the field 

However, there is a much more efficient mechanism for 
producing topological defects. It is due to the facts that 
the field has small fluctuations at the end of inflationary 
stage and that the potential of the field has two minima 
in the radial direction (After inflation, the phase of the 
field $ is almost fixed and it is sufficient to consider only 
the dynamics of the radial direction, and the field has 
small fluctuations beyond the horizon produced at the 
inflationary epoch. See below). In this mechanism, the 
radial part of $ oscillates and settles down to one of two 
minima (±77) of the potential, and even if there are not so 
large fluctuations the final values of $ may be different in 
the different regions in the universe, which leads to the 
formation of topological defects. It is this mechanism 
that we mainly considered in this paper. We search for 
the parameter space where the topological defects are 
produced. 

For a real scalar field, we find that the topologi- 
cal defects are not formed for the GUT scale models 
(i] ~ I0 15 — 10 17 GeV), even if the effective mass square 
of the field m| eff becomes positive. The crucial point is 
that nonthermal fluctuations cannot be as large as the 
amplitude of the coherent mode of the field during pre- 
heating. In fact they are smaller by two orders of mag- 
nitude: (5(f) 2 ) ~ O.O10 2 . Furthermore, even though the 
amplitude of fluctuations becomes the same order as that 
of homogeneous mode due to the effect of the negative 
curvature of the effective potential, the field is dragged 
by its classical motion (homogeneous part) so that the 
field in the entire universe settles down to one minimum 
of its potential. Therefore, in the case of GUT topologi- 
cal defects, the field dynamics (the final value of the field 
$) is determined only by the initial value of its coherent 
mode and the nonthermal fluctuations do not affect deci- 
sively. On the other hand, however, if we are allowed to 
take a breaking scale much smaller than the GUT scale, 
there is formation of topological defects because of the 
initial fluctuations due to inflation. 

For a complex scalar field, fluctuations in the phase 
direction grow much more rapidly and the dynamics of 



homogeneous mode is affected in some extent. This is be- 
cause the field $ does not feel any potential in the phase 
direction. As a result, the amplitude of the fluctuations 
are indeed grows as large as that of the homogeneous 
mode. However, we find that the topological defects are 
not formed in the wide region of the parameter space for 
the GUT scale models. Therefore, we can conclude that 
the defect formation due to the parametric resonance is 
very unlikely, even if the effective mass square of the field 
$ is positive and the amplitude of fluctuations becomes 
large. We also find that topological defects are formed if 
the breaking scale is lower than 10 14 GeV such as axion 
models Q because there are enough time for fluctua- 
tions to affect the dynamics of homogeneous part due to, 
say, the effects of the narrow resonance. 

In Sec. II, we present our models and formalism which 
are used in this paper. We study the real scalar field for 
defects production and find that the GUT defects cannot 
be formed in Sec. III. Section IV is devoted to the study 
of the complex field. Finally, we make our conclusions 
and discussions in Sec. V. 



II. MODEL AND FORMALISM 

We consider the dynamics of a scalar field $. Since 
this field is responsible for the formation of topological 
defects, one can write its Lagrangian as 

c = d^d^*-^m 2 - v 2 ) 2 , (i) 

where g is self-coupling constant and 77 symmetry- 
breaking scale. Here the scalar field $ is real or complex, 
and the produced topological defects are domain walls 
or cosmic strings, respectively. We assume that g is very 
small, otherwise the field $ will reach the minimum of its 
potential very soon during inflation and no topological 
defects are formed. (Here we suppose that the inflaton 
is another field. See below for more discussions.) The 
equation of motion is 

a M ^$ + 3i/$ + 3(|$| 2 -?7 2 )$ = 0, (2) 

where H(= a/a) is the Hubble constant. When we study 
the parametric resonant decay of the complex field $, it is 
more convenient to use two real field X and Y which are 
defined by $ = X + iY with initial conditions ' X(Q) = 
|$(0)| and Y"(0) = 0. We take the initial time at the end 
of inflation. Fluctuations of X and Y can be represented 
in terms of the annihilation and creation operators 

SX(t,x) = — ^— f d 3 k(a x (k)SX k (t)e- rkS 



'In this section, formulas are written assuming that $ is a 
complex field, since the formulas for real <£> is trivial. 
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+al(k)5X* k (t)e ik % (3) 



5Y(t, x) 



(2tt) 3 / 2 



d 3 k(a y (k)$Y k (t)e- ik * 



+aUk)6Y k *(t)e ik *), (4) 



where ai{k) and aj(fc) satisfy commutation relations: 
[eti(fc), aj(fc')] = <5(fc — fc') with i — x,y. If we decompose 
these fields into a homogeneous mode and fluctuations, 
Eq.(0) can be rewritten as 



X + 3HX + g{X 2 -rf)X 

= -3g((5X) 2 )X-g((5Y) 2 )X, 

r k 2 

5X k + 3H6X k + — -grf+ZgX 2 



(5) 



SX k 



= -3g((SX) 2 )5X k - g((SY) 2 )6X k , (6) 



6Y k + 3HSY k 



k 2 

-o - gv 2 + gx 2 



SY k 



mX) 2 )5Y k -Zg{{5Y) 2 )5Y k) 



where (• • •) denotes average over space. The right-hand 
side of these equations represent back reactions due to 
fluctuations and we have used the mean field approxima- 
tions here (e.g., SX 3 ~ 3((5X) 2 )6X ,...). This approx- 
imation corresponds to neglecting rescattering of parti- 
cles , but we can justify it in the following way. Since 
the typical energy of produced particles is estimated as 
E ~ ff 1 / 2 !*! |§, the number density is n ~ g\§\ 4 /E ~ 
^ 1 / 2 |$|3_ Q n the other hand, the cross section of scatter- 
ing is approximately given by 



E 2 
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Then the scattering rate becomes 
r - no ~ 3 3/2 |$| 



(8) 



(9) 



What we have to compare with this is the hubble con- 
stant H ~ X 1 / 2 ^ 2 /M p : 



7 1/2 I$|M„ 



*(*) 



(10) 



where we use g 

i/2|$| _ A i/20 ( see Eq.@ below). 
Therefore, T <C H for g(a(t)/a(0j) <C 1, which means 
that the rescattering effect may be neglected before the 
time when a(t)/a(0) <C g^ 1 . 
Rescaling t,X, and Y as 



"'•Not all the effect of scatterings are neglected in this ap- 
proximations. Actually, forward scatterings are taken into 
account. 



a(r)dr = y/gX(0)a(0)dt, 
Xcl(t) 



V 



X(0)a(0)' J 
Ya(r) 



(11) 
(12) 

(13) 



X(0)o(0)' 

we obtain the rcscalcd dimcnsionless equations 

x" - rfa 2 x + x 3 + 3(5x 2 )x + (6y 2 )x = 0, (14) 
8x1 + [k 2 - fa 2 + ix 2 + 3(Sx 2 ) + (Sy 2 )]Sx k = 0, (15) 
5y'i + [k 2 - f, 2 a 2 + x 2 + (5x 2 ) + 3(Sy 2 )]Sy k = 0, (16) 

where we define k = k/y/gX(0), fj = rj/X(0), and a(0) = 
1. The prime denotes differentiation with respect to r 
and we assume that the universe is radiation dominated, 
which is a reasonable approximation. 

If the universe is radiation dominated, the Hubble pa- 
rameter evolves as 



H(t) = H(0) 



o(0) 
a(r) 



(17) 



Since the potential energy of the inflaton field dominates 
the universe during inflation, the Hubble parameter at 
the end of inflation is 



H(0) 



2ir \/A 



2 (o), 



(18) 



where we assume that the inflaton potential is V(<f>) — 
(A/4)</> 4 (A = 1CP 13 ). First let us consider the case that 
the field $ which we are concerned with is the inflaton, 
i.e., $ = 4> (or g = A). From Eqs.@ and @ 



H(0) 



«(0) 
a(r) 



H(t) = V^(0)a(0)4W 
a 2 {T) 



Integrating this equation and use Eq.(|l8|), we get 

^ (f)(0) 



a(r) 



3 M z 



-o(0)T + a(0). 



(19) 



(20) 



As the inflationary epoch ends at the time when <p = 
Mp/y/Sir 0, and a(0) = 1 is imposed, the evolution of 
the scale factor becomes 



1. 



(21) 



Next we consider the case that the inflaton field is differ- 
ent from <E> field. Then the Hubble parameter becomes 



H(r) = ^X(0)a(0)^=H(0) 

Integrating this equation, we obtain 

, , 0(0) y/\<f>(0) 

a r = ^—t - 

V ' 3 M p ^gX(0) 



o(0) 
a(r) 



(22) 



(23) 



3 



where a(0) = 1 is assumed again. Since the classical 
(homogeneous) mode X slowly evolves as 



X 



1/2 



in the inflationary epoch, Eq.(23) becomes identical to 
Eq.(21) using 0(0) = M p /\/3tt again. Thus we can use 
Eq.(23) for a general scalar field $. 

Now let us consider the initial conditions. The ini- 
tial condition for the homogeneous mode is x(0) = 1 
from its definition. The value of the field $ is initially 
almost constant in the whole universe because of the in- 
flation. However, during inflation, the field has small 
quantum fluctuations, which become seeds for the large 
scale structures of the universe. We take account of these 
small fluctuations in the following way. In the infla- 
tionary epoch, the scalar field fluctuates with the am- 
plitude H/2w(~ 10- 6 (.g/A) 1 / 2 X(0)). These fluctuations 
are stretched beyond the horizon size by inflation and 
become classical. Therefore, one can regard that the ini- 
tial value will be slightly changed as X(0) — > X(0) = 
X{0)(1 + A), where A - 1CT 6 x {g/X) 1/2 , in each region 
of the universe [| and X(0) is considered as the mean 
initial value. Then the initial condition for homogeneous 
mode must be x(0) = 1 + A. 

Next we consider the initial values of fluctuations. 
Hamiltonian is a diagonalized operator for a free field, 
but in the presence of interactions, it is not diagonalized 
in terms of a and a'. It can be diagonalized at any in- 
stant of time by means of Bogoliubov transformations, 
which relate annihilation and creation operators a and 

at t ~ to time dependent annihilation and creation 
operators b(t) and b^(t): 



b x (t) = ax(t)a x + PxttWx, 
b x {t)=f3 x {t)a x +a x (t)a x , 



and a x (t) and (3x it) can be written as 

i J Qx dt 



P. 



e J 



x 



(n x sx k + isx k ), 
(n x sx k - isx k ) 



(25) 
(26) 



(27) 
(28) 



x 



where tt x 



k 2 /a 2 — grj 2 a 2 + 3gX 2 , and the initial con- 
ditions for ax and 0x are 



l<*x(0)| = 1, 



0x(P)=O ) 



(29) 



which means that there is no particles at t = 0. These 
conditions correspond to 



(24) \SX k (0)\ = 



1 



Similarly, initial conditions for SY k are 
1 



, i6X k (0) = n x (0)SX k (0). (30) 



l^(0)| = 



y/2Sl Y (0) 



, iSY k (0) = Q Y (0)5Y k (0), (31) 



where fiy = k 2 /a 2 — grfa 2 + gX 2 . 

Then if we rescale as Eqs.(|IT|)-(|l3|), we get 

1 



|<fe*(Q)| = 
Sx' k (Q) = 

\Sy k (o)\ = 
5y' k (0) = 



y/2Slx(P)X(0y 

Mo) - fi * (0) 



fefe(0), 



(32) 



y/2Q Y [0)X(0) 

M0) 



(33) 



where h(r) = a'(r)/o(T) (from Eqs.@), h(0) = V2/3), 
and the average of the square of the fluctuations can be 
written as PI 



(Sx 2 ) = 
(Sy 2 ) = 



d 3 k 

(27T) 5 

d 3 k 
(2^ 



\5x k \ 2 , 



\Syk\ 



(34) 
(35) 



In the following sections, we integrate Eqs.(|T^)-([l^) 
using relation ( ^T| ) with initial conditions (|32|), ( p3| ) and 
x(0) = 1 + A to see whether topological defects are 
formed or not. We take enough resolution of momentum 
space to calculate Eqs.(|3^) and (35) in all cases which we 
considered. 



III. THE REAL SCALAR FIELD 

A. The classical dynamics of the field: Only with 
cosmic expansion 

First we consider the evolution of the scalar field with- 
out fluctuations after inflation. We expect that the dy- 
namics of the field should be determined entirely clas- 
sically. As mentioned above, however, the scalar field 



s Here the relevant size of the region is between ~ 
H~ l (rj/Alp)' 1 (which corresponds to the horizon at the defect 
formation) and ~ H~ 1 e 60 (which corresponds to the present 
horizon) when |$| ~ 77. 



** Equations (^) and ([35J) should be renormalized pi|| . To 
this end, the vacuum fluctuations are subtracted from these 
expressions when calculated in the numerical calculations. 
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fluctuates with the amplitude H/2ir in the inflationary 
epoch, and the fluctuation will be ~ 10~ 6 x 

{g/X) 1 / 2 at the end of inflation, which leads that the field 
in different space points might fall into the different min- 
ima of its effective potential. Therefore, we should study 
the dependence of the field dynamics (the final value of 
the field) on its initial value. If there are no dependence 
on those initial values whose differences are as large as 
10~ 6 x (g/X) 1 / 2 , we can conclude that the evolution of 
the field is completely determined by classical dynamics. 
Otherwise, we cannot discuss its classical dynamics, to 
say nothing of the effects of nonthermal fluctuations pro- 
duced in the preheating epoch which will be considered 
later. 

From Eq . (|l4|) , the classical equation of motion without 
the effect of fluctuations is written as 



fj 2 a 2 x + x 3 



0, 



(36) 



where the initial condition is x(0) = 1 + A. 

Figure [I] is an example for the time evolution of the 
field with A = (Here we take g = A = 10~ 13 which 
corresponds to the case that the field $ is the inflaton, 
and rj = 10 15 GeV). We see that the amplitude of the field 
x grows larger than the initial value. This is because we 
rescaled the field as Eq.([ll]) and the amplitude grows as 
yfi +i) 2 a 2 {T). 

Figure || shows the numerical result of the critical value 
of A with 77 = 10 12 - 10 17 GeV for g = A = 1(T 13 . If 
A > A cr it, the initial fluctuations do affect the dynamics 
of the field (the final value of the field), which results 
in the formation of topological defects. On the other 
hand, if A < A cr it 5 the fluctuation of the field never 
affects the final value of the field. Aj n f (the amplitude 
of fluctuations at the end of inflation) is less than A cr j t 
for i] > 1 13 GcV. Therefore, the GUT scale topological 
defects are not produced after inflation if we neglect the 
fluctuations of x. 

Let us consider the evolution of the scalar field $ for 
g > A (i.e., $ is not an inflaton). As g grows large, A inf 
grows as Ai„f oc g 1 ^ 2 . On the other hand, the potential of 
the field can be rewritten as V — jiig 1 ^ 2 ^) 2 — (ff 1 ^ 2 ^) 2 ] 2 - 
Since we have relation, g 1/2 $(0) ~ A 1 / 2 0(O) (cf. eq.(f|)), 
the evolution of the field is equivalent to that for g = 
10 -13 when 77 is larger by the factor {g/X) 1 / 2 . Therefore, 
from our calculations (Fig. 0), we find that no topological 
defect is formed for 77 > 10 GeV, which is independent 
of g. 

We can also obtain this result analytically by com- 
paring the initial fluctuation |<5$/$|(~ 1CT 6 x (g/X) 1/2 ) 
with the change of the amplitude AA of oscillating <!> 
due to cosmic expansion in one oscillation time T. Since 
4> ~ — at the critical epoch (when A ~ 77), A A is 
given by 



AA 



HT 



A 1 / 2 ^ 2 



M, 



(!) 



1/2 



where we use eq.(p4|). If AA/A is larger than |(5<f>/<f>|, 
the initial fluctuations do not affect the dynamics of the 
classical evolution of the field <E>. Since A; n f = ~ 
1CP 6 x (g/A) 1 / 2 , we can obtain the condition that the field 
4> settles down to the definite minimum of the potential 
in the entire universe: 



V 



£10" 



(38) 



This is the same as what we obtained from our numeri- 
cal calculations. Notice that this condition is, of course, 
independent of the coupling constant g. 



B. Preheating stage: Effect of nonthermal 
fluctuations 

In order to study the production of nonthermal fluc- 
tuations and find whether the formation of topological 
defects occurs or not, we must investigate how fluctu- 
ations can be created. In fact, they are produced due 
to both parametric resonance during preheating and the 
negative curvature of the effective potential. To inves- 
tigate how large the non-thermal fluctuations grow, we 
consider these two effects separately. First, we study the 
parametric resonance in this subsection. The effect of the 
negative curvature of the potential will be considered in 
the next subsection. 

We integrate those equations which the breaking terms 
are omitted from Eqs. ( |l4| ) and (15), i.e., 77 = 0. For 
the early stage (during preheating), neglecting the break- 
ing terms is quite a reasonable approximation, since the 
breaking scale is much smaller than the initial amplitude 
of the field <f>. The result is shown in Fig. |. The solid line 
denotes the envelope of x 2 and the dotted line 5x 2 : fluc- 
tuations of x, and these lines are extrapolated to larger 
a. Note that the decaying power of homogeneous mode is 
~ —0.26, which is smaller than the value —1/3 estimated 
in the Ref. |l8| . [J This is because the actual dynamics 
is much more complicated since (fa 2 ) is also oscillating 
for example. (Our value —0.26 is much smaller than the 
value ~ —2/3 in Ref. I n this reference, however, the 
calculation has been done for only a short time such as 



A r 



t'The value —1/3 can be estimated with assumptions that 
(fa 2 ) grows large and dominates the effective mass square of 
the field and becomes almost constant. Using this naive es- 
timation, we get the same decaying power for the both cases 
that the evolution of homogeneous mode is written in terms of 
the elliptic function and the trigonometric function. However, 
the dynamics of fluctuation obeys Lame equation, not Math- 
ieu equation, so that the instability band is narrower and the 
growth rate is much smaller in the case of Lame equation. 
Therefore, the decaying power is smaller. 



(37) 
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t ~ 800, so that the decaying power in the narrow reso- 
nance stage cannot be seen from such a short time evo- 
lution.) The time that the field falls into one minimum 
of its potential is estimated as a* ~ (77/ 'X(O))" 1 . For the 
GUT scale (77 - 10 16 GeV), a* - 400. At this time the 
amplitude of fluctuations is much less than that of homo- 
geneous mode. Therefore, the fluctuations is too small 
to produce the GUT defects. On the other hand, we can 
see that the crossing takes place at the time a c ~ 5 x 10 5 . 
For a > a c , the dynamics of the field <£> may be fully af- 
fected by its fluctuations since the field fluctuates more 
than O(l), which means the topological defects may be 
developed in the same sense as pointed out in Ref. p8[ . 
Therefore, IF 77 < 10 13 GeV, it is expected that the full 
symmetry restoration occurs and leads to the formation 
of topological defects due to the later spontaneous sym- 
metry breaking (if there is any breaking term in the equa- 
tion) . 

C. Postpreheating: Effects of the negative curvature 
of the potential 

Now let us take into account the effects of the nega- 
tive curvature of the potential. To this end, we integrate 
equations 

x" - ifa 2 x + x 3 + 3(Sx 2 )x = 0, (39) 
5x'l + [k 2 - 77V + 3x 2 + 3(Sx 2 )]Sx k = 0, (40) 

and see how large the amplitude of the fluctuations grows 
comparing with that of the homogeneous mode. The 
time evolution of the homogeneous mode and the fluctu- 
ations are shown in Figs. [| and ^, respectively. From 
these figures we see that the preheating stage ends at 
a ~ 250. What we must pay attention to is the fact that 
the amplitude of nonthermal fluctuations during preheat- 
ing is two orders smaller than that of homogeneous mode, 
which is contrary to the estimation of Ref. jl8| (It can 
also be seen from Fig. ||, and is consistent with Ref. ) . 
Moreover, the amplitude grows due to the global instabil- 
ity of negative mass (breaking term), but the maximum 
amplitude is at most a few factor less than that of the 
homogeneous mode (at a ~ 3000). This is the crucial 
point. We can say in other words as follows: although 
the effective mass square of the field becomes positive 
during preheating (see Fig. ^), the amplitude of the fluc- 
tuations is smaller than that of homogeneous one so that 
no full symmetry restoration takes place, which results 
in no development of GUT topological defects (domain 
walls). Notice that the typical momentum of the fluc- 
tuations which contributes to fluctuations (Sx 2 ) is larger 
than H ~ \f\cf) ~ v / g| ( l > | (see Fig. 0). This means that it 
is a good approximation to take the conformal vacuum 
for initial value of fluctuations. 



D. Dynamics of the field with the initial fluctuations 

As we mentioned in Sec. II, there is much more efficient 
mechanism for producing topological defects. It is due to 
the facts that the field has small fluctuations at the end of 
inflationary stage and that the potential of the field has 
two minima in the radial direction. In this mechanism, 
the radial part of $ oscillates and settles down to one 
of two minima (±77) of the potential, and even if there 
are not so large fluctuations the final values of <i> may be 
different in the different regions in the universe, which 
leads to the formation of topological defects. In order to 
study this effect is the main investigation of this paper, 
and we show it in this subsection. 

We consider the dynamics of the field taking into ac- 
count of the initial fluctuations created during inflation 
in the following way: As we discussed in Sec. II, we nu- 
merically integrate Eqs. ( |39"| ) and (|i(]) with initial con- 
ditions ([52j) for fluctuations and x(0) = 1 + A for a ho- 
mogeneous mode, where A denotes the amplitude of the 
initial fluctuations, and see how much the dynamics of 
the homogeneous mode is affected by those nonthermal 
fluctuations produced both during preheating due to the 
parametric resonance and due to the negative curvature 
of the potential. 

We can see little effects of the nonthermal fluctuations, 
since the difference between Figs. || and ^ cannot be 
seen. Therefore, we conclude that topological defects are 
not formed for 77 > 10 13 GeV and any value of g, which 
means that the GUT defects are not be produced during 
reheating after inflation. This is the same conclusion as 
the case without considering the nonthermal fluctuations. 

IV. THE COMPLEX SCALAR FIELD 

A. Preheating stage: Effect of nonthermal 
fluctuations 

Let us turn our attention to the dynamics of a com- 
plex field <& = X + iY. As mentioned previously, there 
are three kind of fluctuations: the nonthermal fluctua- 
tions created due to parametric resonance, the one due 
to the effect of the negative curvature of the potential, 
and the initial fluctuations produced during inflation. In 
this subsection, we pay attention to the first one. To this 
end, omitting the breaking term, we have integrated the 
following rescaled equations: 

x" + x 3 + i(5x 2 )x + (Sy 2 )x = 0, (41) 
6x k + [k 2 + 3s 2 + S(Sx 2 ) + (Sy 2 )]Sx k = 0, (42) 
5y' k ' + [k 2 +x 2 + (Sx 2 ) + 3(Sy 2 )]Sy k = 0. (43) 

The result is shown in Fig. g. The solid line denotes the 
envelope of (p 2 = x 2 and the dotted line fluctuations, i.e., 
(Sip 2 ) = (Sx 2 ) + (Sy 2 ), which are extrapolated to larger a. 
At a*(~ (7j/X(0)) _1 ) ~ 400, there are not large enough 



G 



fluctuations to produce the GUT defects. We see that 
the crossing time is a c ~ 10 5 . This means that topologi- 
cal defects may be formed for the models with breaking 
scale r\ < 3 x 10 13 GeV. Note that the amplitude of fluc- 
tuations is larger than that of the real scalar field. This is 
because the phase fluctuation 59 (Goldstone mode) ap- 
proximately written as Sy/ x feels no potential so that 
it can grow much faster po| , po|| . Therefore, the crossing 
time is earlier than the real scalar field case even if the 
decaying power (~ —0.185) is smaller due to large back 
reaction from the phase fluctuation. 



B. Postpreheating: Effects of the negative curvature 
of the potential 

Now we must include the breaking term in the equa- 
tions in order to see the effect of the negative curvature 
of the potential on producing fluctuations. We integrate 
Eqs.(^) - ([!(]) numerically. The dynamical evolution 
of the homogeneous mode with taking account of the 
effects of non-thermal fluctuations is shown in Fig. [r]. 
The field decays abruptly at a ~ 50. This time is ear- 
lier than that for the case of the real scalar field. The 
reason can be seen in Fig. Il2| (Figs. 
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and |12| are the 
evolutions of radial and phase fluctuations, respectively, 
when A = 0). It is reasonable that the amplitude of Sy 
grows much faster than that of Sx (compare Fig. [ll] with 
Fig. ^) , since the fluctuation Sy has an approximate re- 
lation S9 ~ Sy/x to the phase fluctuation which can be 
identified with a Goldstone mode who feels no potential, 
as previously mentioned. Therefore, the energy stored 
in the homogeneous mode of the field was quickly trans- 
ferred into this Goldstone mode first. This is also seen in 
the momentum distribution of Sx and Sy (Figs. [r| and 

We see, however, that nonthermal fluctuations are not 
so much produced in the preheating epoch. At the end 
of preheating, the amplitudes of nonthermal fluctuations 
are ((Sx) 2 ) ~ 0.01a; 2 , ((Sy) 2 ) ~ O.lx 2 , which are smaller 
than those estimated (~ x 2 ) in the Rcf. p8| . There- 
fore, topological defects cannot be directly formed due 
to the effects of nonthermal fluctuations produced, even 
if the effective mass square of the field is positive during 
preheating (see Fig. [l5|). Even though the amplitude of 
fluctuations becomes as large as that of the homogeneous 
mode at a ~ 3000 due to the effect of the negative curva- 
ture of the potential, it is still not enough condition for 
the formation of topological defects (It will be verified in 
the next subsection). 



C. Dynamics of the field with the initial fluctuations 

As mentioned many times above, there are fluctuations 
produced at the inflationary epoch. Therefore, the ini- 
tial value of the field $ is slightly different in each region 



which contains many horizons, and we take the initial 
condition as x(Q) = 1 + A. In order to consider the ef- 
fects of nonthermal fluctuations to the dynamics of the 
homogeneous mode in the presence of this initial fluctua- 
tions, we integrate Eqs.([l4) - ( |l6| ) with initial conditions 
@ and (||) for fluctuations and x(0) = 1 + A for the 
homogeneous mode. 

Figure [r] is the result for g = 10~ 13 . The evolution 
of the homogeneous mode (the final value of the field) 
is affected if the breaking scale ij < 10 14 GeV. Compar- 
ing with the results without including fluctuations (the 
dashed line), the fluctuations considerably affect the dy- 
namics of the homogeneous field (the final value of the 
field), and it is different from the real scalar field case. 
This is because the large phase fluctuations are produced. 
There are some unstable region about r\ ~ 10 16 GeV. For 
r\ ~ 10 16 GeV, the critical time when the field settles 
down to one of the minimum of the potential is almost 
the same time when the radial fluctuations become sig- 
nificant. We think that this is the reason for the unstable 
behavior of A around r\ — 10 16 GeV. 

For g > A, we can study the dynamics of $ by rescaling 
Ai„f and rj from the result for g = A = 10~ 13 . In Fig. [l7] 
we show the region of the defects formation in the g — rj 
plane. From this figure it is seen that the topological 
defects are not formed for almost all the region of the 
GUT scale (rj = 10 15 - 10 17 GeV). 

Nonthermal fluctuations do actually affect the evolu- 
tion when comparing with the results for the classical 
evolution without fluctuations studied in the previous 
section. Owing to the effects of nonthermal fluctuations, 
the region where topological defects are produced are 
larger than that for the pure classical one. Note that the 
GUT symmetry is not fully restored during preheating in 
spite of the fact that the effective mass square of the field 
becomes positive before the field settles down to the min- 
imum of its effective potential. This is because that the 
amplitude of fluctuations do not highly overcome that 
of the homogeneous mode, and the fluctuations cannot 
disturb the whole dynamics of the field even if its am- 
plitude becomes as large as that of homogeneous mode. 
The authors in Ref. |l8j assert that the criterion of sym- 
metry restoration is the sign of the effective mass square 
of the field. As we mentioned above, this condition is 
not enough to judge the symmetry restoration. They 
consider that the amplitude of fluctuations is as large as 
that of homogeneous mode at the end of preheating and 
the former always overcome the latter due to the nar- 
row resonance effects in the following stage. The very 
difference from their conclusion is that the fluctuations 
of the filed $ are not created so much and hence there 
needs much longer period of narrow resonance stage for 
the amplitude of fluctuations to be much larger than that 
of homogeneous one. 

Another viewpoint is insisted by the authors of Ref. 
p2fl : the symmetry is always restored if there is a very 
large energy density in the initial state enough for the 
field to occupy both vacua under the dynamical evolu- 
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tion. However, we see that the topological defects will 
not be formed even if initial energy is large. Therefore, 
the criterion given in Ref. 0] is also not enough condi- 
tion. 



V. CONCLUSIONS AND DISCUSSIONS 

We consider the dynamics of both the real and complex 
scalar field $ with potential (g/2)(|$| 2 — rj 2 ) 2 after infla- 
tion. In the preheating phase, the homogeneous mode of 
the field $ transfers its energy into nonthermal fluctua- 
tions through the parametric resonance, and their ampli- 
tudes are expected much larger than that in the thermal 
equilibrium. These nonthermal fluctuations may cause 
the restoration of symmetry and hence the formation of 
topological defects. In the present work, we consider tree 
types of fluctuations. The first one is produced due to 
parametric resonance effects, the second is due to the ef- 
fect of the negative curvature of the potential, and the 
last is created during inflation. 

To investigate these fluctuations thoroughly, we nu- 
merically integrated the coupled equations for a homoge- 
neous mode and fluctuations using the mean field approx- 
imations for the very flat potential U($) = (g/2)(\&\ 2 — 
rj 2 ) 2 . At the inflationary epoch, the field $ fluctuates 
and its amplitude becomes 8§ = H/2tt. Therefore, at 
the end of inflation, <5$/$ - 1CT 6 x (g/X) 1/2 . This ini- 
tial fluctuation might affect the dynamics of the field (the 
final value of the field, i.e., $ = 77 or —77) in the reheating 
epoch after inflation. On account of the initial fluctua- 
tions, we take the initial value of the homogeneous mode 
such as x(0) — 1 + A, which is equivalent to consider- 
ing each domain of the universe beyond the horizon with 
different initial conditions. 

First we study the evolution of a classical field with 
various initial conditions A when there is no nonthermal 
fluctuations produced by the effect of parametric reso- 
nance, and find that the topological defects are formed 
for 77 < 10 13 GeV. The physical meaning of this is very 
clear. If the amplitude of initial fluctuations is smaller 
than the amplitude change in one oscillation time due 
to cosmic expansion, it has no influence over the field 
dynamics. Since the energy change can be estimated as 
A A/ A ~ H/uj osc ~ Mqut /Mp ~ 10' 4 for the GUT 
models, <C AA/A so that the initial fluctuations 

cannot affect. 

Next we study the dynamics of the real scalar field tak- 
ing into account the effects of nonthermal fluctuations 
produced as a result of both the parametric resonance 
and the effect of the negative curvature of the potential 
in the preheating epoch, and how they affect the dynam- 
ics of the field. We find that fluctuations do not influence 
the dynamics of the homogeneous mode at all, and it is 
entirely the same as the case of the classical evolution 
without fluctuations. Therefore, no GUT topological de- 
fects (77 ~ I0 16 GeV) can be formed even if the amplitude 



of the fluctuations produced during the reheating stage 
grows as large as that of the homogeneous mode. This 
means that the symmetry restoration does not occurred 
in spite of the positiveness of the effective mass square of 
the field $. 

For the complex scalar field, fluctuations do affect the 
evolution of the field, since the phase fluctuation is pro- 
duced much more rapidly than that in the radial direction 
since the field feel no potential in phase direction (Gold- 
stone mode). However, there is no but a tiny region for 
the formation of topological defects for the GUT scale 
in the parameter space, even if fluctuations are much 
produced and the effective mass square of the field be- 
comes positive during preheating. In other words, it is 
not enough condition for the symmetry restoration that 
the field has a positive effective mass square and fluctu- 
ates at most of order 0(1). On the other hand, topo- 
logical defects are formed in the model with much lower 
breaking scale (77 < 10 14 GeV) such as axion models. 
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FIG. 1. Evolution of the scalar field without nonthermal 
fluctuations (only classical evolution) for the initial deviation 
A = and g = 10" 13 . 

FIG. 2. Critical value of A as a function of r) for g = 10~ 13 . 
The stars are from numerical calculations. The dashed line 
shows the analytic estimation of Eq.(^) including some nu- 
merical factors, and the dotted line denotes the amplitude of 
the initial fluctuations produced in the inflationary epoch. If 
the latter line is above the former one, the dynamics of the 
field is considerably affected, which results in the formation 
of defects. 

FIG. 3. Long time evolution of the dynamics of the real 
scalar field with nonthermal fluctuations. The dots and stars 
denote the amplitude square of homogeneous part (a; 2 ) and 
fluctuations ((Sx 2 )), respectively. The solid line stands for the 
envelope of x 2 and is extrapolated to the large a. Similarly, 
the dotted line represents the envelope of (8x 2 ). 

FIG. 4. Evolution of the homogeneous mode of the real 
scalar field with nonthermal fluctuations with initial deviation 
A = 0. We take rj = 10 15 GeV and g = 10" 13 . 

FIG. 5. Evolution of the fluctuation of the real scalar field 
with initial deviation A = 0. 

FIG. 6. Evolution of the effective mass square at <E> = 
for the real scalar field with initial deviation A = 0. 

FIG. 7. Power spectrum of the radial fluctuation 5xt for 
(a) a ~ 1000, (b) a ~ 2000, (c) a ~ 3000 and (d) a ~ 3500. 

FIG. 8. Critical value of A as a function of r\. The stars 
are from numerical calculations. The dashed line shows the 
analytic estimation without fluctuations (see Fig. ^), and the 
dotted line denotes the amplitude of the initial fluctuations 
produced in the inflationary epoch. 

FIG. 9. Long time evolution of the dynamics of the com- 
plex scalar field with nonthermal fluctuations. The dots and 
stars denote the amplitude square of homogeneous part (<f> 2 ) 
and fluctuations ((8cj> 2 )), respectively. The solid line stands 
for the envelope of 4> 2 and is extrapolated to the large a. Sim- 
ilarly, the dotted line represents the envelope of {8(j> 2 }. 

FIG. 10. Evolution of the homogeneous mode of the com- 
plex scalar field with nonthermal fluctuations for initial devi- 
ation A = and g = 10~ 13 . 

FIG. 11. Evolution of the radial fluctuation 8x of the 
complex scalar field with initial deviation A = 0. 

FIG. 12. Evolution of the phase fluctuation 8y of the 
complex scalar field with initial deviation A = 0. 

FIG. 13. Power spectrum of the radial fluctuation 5xk for 
(a) a ~ 1000, (b) a ~ 2000, (c) a ~ 3000 and (d) a ~ 3500. 

FIG. 14. Power spectrum of the fluctuation 8y^ for (a) 
a ~ 1000, (b) a ~ 2000, (c) a ~ 3000 and (d) a ~ 3500. 

FIG. 15. Evolution of the effective mass square at $ = 
for the complex scalar field with initial deviation A = 0. 
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FIG. 16. Critical value of A as a function of rj. The stars 
are from numerical calculations. The dashed line shows the 
analytic estimation without fluctuations (see Fig. ^), and the 
dotted line denotes the amplitude of the initial fluctuations 
produced in the inflationary epoch. 



FIG. 17. Parameter space for formation/ no formation of 
topological defects in the g — r\ plane. NA denotes the region 
where we do not have the result, and OC is the region where 
the initial value of the field is less than the breaking scale: 
|$(0)| < 77. In the OC region, the field already settles down 
to one the minimum of the potential at the end of inflation. 
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